Spin-glass instability of short-range spherical ferromagnet 
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In structurally disordered ferromagnets the weak random dipole-dipole exchange may transform 
the polydomain state into a spin-glass one. To some extent the properties of such phase in disordered 
isotropic ferromagnet can be qualitatively described by the spherical model with the short-range 
ferromagnetic interaction and weak frustrated infinite-range random-bond exchange. This model 
is shown to predict that spin-glass phase substitute the ferromagnetic one at the arbitrary small 
disorder strength and that its thermodynamics has some similarity to that of polydomain state 
along with some significant distinctions. In particular, the longitudinal susceptibility at small fields 
becomes frozen below transition point at a constant value depending on the disorder strength, while 
the third order nonlinear magnetic susceptibilitiy exhibits the temperature oscillations in small field 
near the transition point. The relation of these predictions to the experimental data for some 
disordered isotropic ferromagnets is discussed. 

PACS numbers: 64.60.Cn, 05.70.Jk, 64.60.Fr 

The spherical model with short-range exchange shares the basic qualitative features with real isotropic ferromagnets. 
It has phase transition only in space dimensions d > 2 and its coercive field is strictly zero. This is because it is 
equivalent to the isotropic n-component model in the limit n — > oo So the scalar magnetization of short-range 
spherical model corresponds to the magnetization module of isotropic ferromagnets and this makes this model very 
useful for the studies of qualitative features of their thermodynamics. Yet real ferromagnets have also the long- 
range dipole-dipole interaction. Being a weak relativistic effect it nevertheless determines crucially the nature of 
ferromagnetic transition which usually results in the appearance of inhomogeneous polydomain state. It shows up in 
the freezing of longitudinal magnetic susceptibility at the value x — (47tk) _1 below T c at fields H < Ai:k,M s , k is the 
, depolarizing coefficient along the field direction, M s is the spontaneous magnetization @, Q . It is rather natural to 
suppose that when some non-magnetic disorder such as structural defects or non-magnetic impurities is present in a 
O ! crystal the polydomain state may transform into the spin-glass one 

To describe the qualitative features of such spin-glass state in random isotropic ferromagnets we may turn to 
the spherical model with weak long-range frustrated disorder imitating the random dipole-dipole exchange in the 
structurally disordered media. The influence of such (infinite-range) disorder on the thermodynamics of the mean- 
k*" ' field spherical ferromagnet was studied in Ref. [H In this model the spin-glass phase instead of ferromagnetic one 
do appear when disorder becomes sufficiently strong. Here we consider more realistic short-range spherical model of 
ferromagnet with the same infinite-range frustrated random exchange. We find that contrary to the mean-field model 
in the short-range one the spin-glass substitutes the ferromagnetic phase at arbitrary weak random exchange. We 
• ■ also show that the magnetic properties of this spin-glass phase in small magnetic fields have some similarity to those 
of polydomain ferromagnetic state along with some significant distinctions. 
. The Hamiltonian of the spherical model has the form 
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Here S is TV-component vector subjected to the constraint S 2 = N, H is the external field and Jj j is the matrix of 



-,3 

^ . exchange integrals. 

Partition sum of the model can be represented as 
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(3 = 1/T is the inverse temperature. The parameter a in the integral over A can be arbitrary provided it obeys the 
condition a > /3J max , J max being the largest eigenvalue of J. Thus Eqs. |l|[3]) are valid for any J with the spectrum 
limited from above. For the equilibrium thermodynamic potential F we have from Eq.([T]) at N — > oo 



F = mm F(X)=F |A (yJJ . (4) 

Here Ao ^J^ is the value which provides the minimum of F (A). It obeys the equation of state 

dF(Xo) 
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Solving Eq. ([5]) and substituting the Ao ( J^j found into Eq. ([3]) we get the equilibrium potential F and can then find 
all thermodynamic variables of the system. In particular, we get for the average local spins 

<S) T = -iv|| = /3G(Ao)H (6) 

When J is a random matrix we should average F over it. It can be easily done if we assume Ao ^J^ to be the 
self-averaging quantity. Then while averaging of Eqs. (4, 5) we can just substitute Ao ^j) by its average value 
X=(Xo(f)^y Thus we get from Eqs. (@ gj) 

-2pF=-2p(F(\))j = \- J dep(e)ln{X-e)+N- 1 (3 2 H^G(X)^ ^H, (7) 
where p (e) is the average spectral density of the matrix (3 J, 

p(e) = — lim ImTr (6(e-iS)) . (8) 
From Eqs. ^ [31 [5]) we get the equation for A, 

D (A) + Q (A) - 1, (9) 

D(X)=N- 1 Tr(d(X)) } , (10) 

Q(A) =N- 1 {i 2 K(G 2 (\))K = N- 1 ((S) 2 T ) j (11) 

The last equality in Eq. (fTl]) follows from Eq. © . It shows that Q (A) is the Edwards- Anderson spin-glass order 
parameter. 

Here we consider the Gaussian disorder for the exchange integrals with the mean 

(Ji,j) = J (r< - Tj) 

and the deviation 

A 2 

We assume J(r, — r^) to describe the short range ferromagnetic interactions so its Fourier transform J(k) have a 
maximum at k = and near it 

J(k) « J-Ak 2 . 

Then on a three-dimensional lattice the spectral density of (3 J (k), 
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would have the square-root behavior at the upper edge of the spectrum which describes the most relevant long-range 
ferromagnetic fluctuations, 



So we choose 
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Here 9 is the Heaviside step function. This po (e) correctly behaves at the upper edge and makes further calculations 
quite easy. The explicit form of J (k) appears to be irrelevant for the homogeneous external field we consider below 
and all thermodynamics is determined solely by po (e). 

Now we can find ^G (A)^ for such random ensemble where the weak infinite-range random exchange fluctuations 

of arbitrary sign coexist with non-random short-range ferromagnetic interactions. Expanding G (A) in the power 
series of J[ ^ and averaging this expansion with the Gaussian distribution we find in the large N limit the following 

expression for the Fourier transform of ^G (A)y , 

G _1 (A,k) = A-/3 2 A 2 L>(A) -/3J(k). (13) 
Then for D (A) (fT0|) we have the equation 
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The solution to this equation is 
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Eqs. H3HH1) define G (A, k). From these equations we can also find the Fourier transform of (G 2 (A) 



(14) 



G 2 (A, k) = -jLG (A, k) = G 2 (A, k) [1 - (i 2 A 2 D< (A)] 



From dHl EE HIl we also get 
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Thus we have all that is needed to obtain the explicit expressions for the average thermodynamic potential ([7]) and 
the equation of state ©■ Further we consider the homogeneous external field, Hi — H,i — 1, , N. It is convenient to 
introduce the new variable z, < z < 1, instead of A, 



t PJ f -1 \ 



Then we have from Eqs. (9-11, 13-15, 17) the equation which defines z, 

h 2 z 2 (1 + cz) = (1 - tz) (1 - z 2 ) (1 - czf 



(17) 
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h = H/Tg, T g = ^JJ 2 /A + A 2 ,t = T/Tg. 
From Eqs. ([H HH [TH [17]) we get the average potential, 

- 2F/T g = tint + z + z~ 
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It can be easily checked that Eq. (fTH|) is equivalent to the equation §§- = and that the solution of it provides the 
minimum of potential in the interval < z < 1. Other thermodynamic parameters can be also expressed via z. Thus 
averaging Eq. (|6|) over random exchange we get the average magnetization 



M = (3HG (A, k = 0) 



z h 



\2 : 
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while from |9j QUI flT|) we get for the equilibrium value of the Edwards- Anderson order parameter 

Q = l-tz. 

Also from Eq. (fTi?)) we obtain the entropy 
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and the heat capacity 
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Thus Eqs. (|20j|23[) supplied with the solution to Eq. (|18[) for z = z (t, h, c) give full description of the thermodynamics 
of the model. Here we should note that parameter c defined in Eq. (| 14[) determine the relative strength of the 
short-range ferromagnetic bonds. It varies in the interval < c < 1 and c = 1 corresponds to the pure short-range 
ferromagnet while at c = only random infinite-range glassy exchange is present in the system. So at c = 1 we 
have ordinary ferromagnetic transition at t = 1, h = with anomalies usual to the pure spherical model. In this case 
Q = M 2 . Yet at all c < 1 this transition is destroyed and instead the transition into the spin-glass phase takes place 
at t = l,h= 0. 

Indeed, when h — > then z — ► 1/t for t > 1 and z — *■ 1 for t < 1. So at all c < 1 and ft. = M is zero, but 
spontaneous Q appears at t < 1, Q = 1 — t. This is in sharp contrast with the model where instead of short-range 
J(k) the infinite-range mean-field ferromagnetic interaction of the form J(k) = <5k.o^/A is introduced Then 
spin-glass transition substitutes the ferromagnetic one only at c 2 < 1/5. 

As the short-range spherical ferromagnet correctly reproduces the qualitative features of real isotropic ferromagnets 
we may suppose tentatively that the destruction of magnetic order by the infinitesimal glassy long-range random-bond 
disorder can also take place in real isotropic magnets with structural imperfections. In such a case the present model 
can reveal the qualitative features accompanying this phenomenon in some amorphous ferromagnets, ferromagnetic 
alloys and even in the nominally pure ferromagnetic crystals. 

The dependence of the spin-glass transition temperature T g (18) on the relative strength (A/ J) of frustrated disorder 
is shown in Fig. [TJ 
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FIG. 1: (color online) The dependence of T g on the relative disorder strength A/J . 
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FIG. 2: (color online) Field dependence of magnetization for c = 0.9 at various temperatures, t = 3 (solid line), t = 1 (dashed 
line), t — 0.1 (dotted line). 



So the zero-field thermal properties of the model do not depend on c, while the magnetic susceptibility % is essentially 
defined by it. Fig. [2] shows the field dependence of magnetization at various temperatures. Note the steep rise of 
M at low fields below T g . Here the slope of M(h) in small fields is limited by the value of zero-field susceptibility, 
(1 — c)~ 2 , while in the polydomain ferromagnet it is limited by \ — (47tk) _1 . We can easily find M(t, h) from Eqs. 
d [201) for small fields h 2 < t (1 - c) 3 , 



M = h 



(1 - c) 2 + (1 - c 2 ) It + Vr 2 + bh 2 



t - 1 
2t 



b = 



2t (i - c y 



(24) 



Fig. [3] presents this M(t) and \{t) m small fields for c = 0.9. They are rather similar to those of pure ferromagnet 
undergoing the transition into polydomain state albeit with the disorder-dependent saturation values. 




FIG. 3: (color online) Temperature dependencies of M (a) and \ 3 ) f° r c = 0-9 in small fields, h — 5 x 10 4 (dotted lines), 
3 x 10" 4 (dashed lines), 10" 4 (solid lines). 

Yet more spectacular anomalies are exhibited by the nonlinear magnetic susceptibilities of the model. They are 
known to diverge at spin-glass transition in zero field in various mean-field spin-glass models 6] including the spherical 
one In the last case these divergences result from the specific non-analyticity of M(t, h, c) at t = 1, h = in Eq. 
(|24|) which also give rise to temperature and field oscillations of nonlinear susceptibilities near the transition. Near 
transition point at c ^ 1 and for |r| 1 — c we get from Eq. (|24p two first nonlinear magnetic susceptibilities, 

d 2 M _ 2b 2 h (3t 2 + 2bh 2 ) d 3 M _ 66 2 r 4 

X2 = ~ {T 2 +bh 2f2 ^3 = - — - bh2)5/2 ( 5) 
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They exhibit highly anisotropic behavior near the singular point r = 0, h = 0. In the polar coordinates defined as 

r = \J t 2 + bh 2 cp = tan -1 f J 



we have 



X2 



6b 2 

26 3/2 sin ip (2 + cos 2 , X3 = cos 4 



(26) 



Thus at <p = 0(h = 0) X2 = 0, X3 = 6& 2 |r| 1 , while at (p = n/2(r = 0) \2 = 2b 3 / 2 sign(h), \z = 0. The behavior of 
Xi and X3 near the singular point r = 0, h = is shown in Fig. [4] 
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FIG. 4: (color online) Nonlinear susceptibilities x^/b 3 ^ 2 (a) and x^/b 2 (b) near the singular point r = 0, /i = 0. 




FIG. 5: (color online) (a) - field dependencies of x^/b 3 ^ 2 at r = 0.01 (solid line), 0.03 (dashed line), 0.05 (dotted line); (b) - 
temperature dependencies of Xzjb 2 at small fields h = (solid line), 0.001 (dashed line), 0.002 (dotted line). 

These complex anomalies result in specific field dependence of xi and temperature oscillations of X3 as seen in Fig. 

El 

The behavior of nonlinear susceptibilities similar to that of Fig. [5] is observed in isotropic ferromagnet 
Ndo.75Bao.25Mn03 @ and in the polycrystalline samples of RuSr2GdCu20s^- In the toroidal polycrystallinc 
samples of Lao^eBao^MnO^ with demagnetization factor k ~ the plateau in x(T) same as in Fig. [3]Jb) is found 
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manifesting the transition into the glass state |10j. There are many other examples of such step-like behavior of x(T) 
in disordered isotropic magnets, see, for example, Refs. [TTJ] , [l^. But it is often impossible to check the relation 
X = (47tk) _1 below T c to distinguish between the polydomain and the spin-glass states as some experimental papers 
lacks the values of k calculated from the sample shape. It is quite possible that such check will show that many 
allegedly polydomain ferromagnets are actually the spin-glasses. 

Yet now it is not clear if the present result on the spin-glass instability of spherical magnet does apply to the 
real dipolar Heisenberg magnets which may have some threshold disorder strength to become the spin-glasses. To 
resolve this issue further theoretical studies of the role of random dipole-dipole interaction in isotropic ferromagnets 
are needed. 

I gratefully acknowledge the useful discussions with V.B. Shirokov and E.D. Gutlianskii. 
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